Abstract. Recent observations of Krucker & Benz (1998) 



give strong support to Parker's hypothesis ( 1988 ) that 
small scale dissipative events make the main contribu- 
tion to quiet Sun coronal heating. They also showed that 
these small scale events are associated not only with 
the magnetic network, but also with the cell interiors 



(Benz & Krucker, 1998| ). Taking into account in addi- 



tion the results of the analysis performed by Priest with 



co-authors (2000) who demonstrated that the heating is 
quasi-homogeneous along the arcs we come to the con- 
. . , elusion that the sources driving these dissipative events 

• are also small scale sources. Typically they are of the or- 
| der of or smaller than the linear scale of the events ob- 

^vq ■ served, that is < 2000 km. To describe statistical proper- 
ties of quiet Sun corona heating by microflares, nanoflares, 
and even smaller events, we consider a cellular automata 
model subject to uniform small scale driving and dissi- 
pation. The model consists of two elements, the magnetic 
C \ field source supposed to be associated with the small scale 
hydrodynamic turbulence convected from the photosphere 
and local dissipation of small scale currents. The dissipa- 
^ ■ tion is assumed to be provided by either anomalous resis- 
tivity, when the current density exceeds a certain thrcsh- 
\l ■ old value, or by the magnetic reconnection. The main 
\ problem considered is how the statistical characteristics 
of dissipated energy flow depend upon characteristics of 
the magnetic field source and on physical mechanism re- 
f^) i sponsible for the magnetic field dissipation. As the thresh- 
old value of current is increased, we observe the transition 
O from Gaussian statistics to power-law type. In addition, 
' | we find that the dissipation provided by reconnection re- 
^ ■ suits in stronger deviations from Gaussian distribution. 

• Key words: solar corona - microflares - heating - cellular 
^ automata model 
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1. Introduction 

Although photospheric convection may supply enough en- 
ergy to heat the coronal plasma, the way how energy is 
supplied and dissipated in t he cor ona is not complet ely un - 
derstoo d (e.g . Priest et al., |2000| ; Einaudi & Velli, |l994| ) 



Parker ( 1988] ) suggested that the solar corona could be 



heated by the dissipation at many small scale tangential 
discontinuities arising spontaneously in the coronal mag- 
netic fields braided and twisted by random photospheric 
footpoint motions. To emphasize these events, he intro- 
duced the notion of nanoflares. This idea stimulated the 
intensive search of observational signatures of microflares 
and nanoflares as well as many theoretical developments 
on the contribution of small scales to energy dissipation 
in the solar corona. The solar flares, most energetic im- 
pulsive phenomena in active regions were observed and 
studied in detail for quite a long time. On a large scale, 
it is well known that the most energetic dissipative events 
occur in the vicinity of tangential discontinuities at the 
borders of the photospheric dipolar regions like bright 
spots. These sources create strong currents resulting in 
reconnection and flare-like events. Such eruptive events 
are related to magnetic energy releases, sudden changes of 
the magnetic field topology, and plasma heating. There- 
fore Parker's idea about solar corona heating by nanoflares 
used quite similar physical concepts on the nature of dis- 
sipative phenomena, but it pointed out the importance of 
small scales. 

The smaller scale phenomena, microflares, were first 
detected in soft X-rays in a balloon experiment by Lin et 
al. ( |1984[ ). The development of new instrumentation al- 
lowed to perform multi-wave satellite (YOKOH, SOHO, 
TRACE) and ground based (VLA) high-resolution obser- 
vations of smaller scale (about thousand kilometers) lower 
energy phenomena. They were observed in active regions 
but also in the quiet regions of the Sun and in coronal 
holes. The small intensity flares observed in active re- 
gions were quite rare and not powerful enough for heat- 
ing (Shimizu et al., 1994). These observations were per- 



formed in the regions of high level of background flux and 
strong stray light. The very same microflares can be more 
easily detected in the quiet corona. Soft X-ray observa- 
tions (Benz et al., 1997 ) and EUV observations (Harrison, 



1997) have revealed enhanced emission and thus intense 



heating above the magnetic network. Similar phenomenon 
that form small X-ray jets at the limb was reported by 
Koutchmy et al. ( 1997 ). It was found that the number of 
observed events increases with the sensitivity. 

Another way to shed some light on the important prob- 
lem concerning the characteristic scales of the major dis- 
sipative events consists in statistical study of different ob- 
served parameters. The major objective of such a study is 
to obtain the Probability Density Function of energy in a 
wide range of energies. 

For regular flares that occur mainly in active regions 
Crosby et al. ( 1993 ) have found that the total energy in 
the flare electrons observed in hard X-ray bremsstrahlung 
has a power law distribution f(E) ~ E a with index a = 
— 1.53± 0.02. But the energy supplied by the flares in the 
active regions is not sufficient for the corona heating. 

The crucial factor that determines the characteristics 
of the heating process and its scales is related to the power 
law index of the PDF. If it is larger than minus two the 
major heating takes place in the intensive large scale dis- 
sipative events. If it less than minus two it occurs in small 
scales. The reason is that to provide the explanation of 
the efficient coronal heating in small scales it is necessary 
to have an excess in the occurrence rate of small scale 
events. An important result that supports the hypothesis 
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of Parker was reported by Krucker & Benz (199S) who 
have found from the Yohkoh / SXT observations that the 
PDF dependence upon the energy has a power law distri- 
bution in the energy range 10 24 -10 26 ergs with the index 
about —2.59. This result was obtained assuming that the 
flaring region has a constant height. The authors not only 
concluded that the weak flaring events rule the heating 
process, but they also noticed that the brightest small 
scale microflares lie above enhanced elements of the mag- 
netic network, and the denser ones exhibit a higher level 
of fluctuations. This means that the corona is continu- 
ously replenished by underlying chromospheric material 
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that has be en heated to coronal temperatures. Mitra & 
Benz ( 2000 ) have discussed the same observations but sup- 
posing that the height varies proportionally to the square 
root of the area and have shown that the index becomes a 
little larger but still is smaller than minus two. This result 
was confirmed by Parnell & Jupp ( 2000 ), who estimated 
the index to be between —2 and —2.1 making use of the 
data of TRACE. However, Aschwanden et al. (2000) using 
also the data of TRACE have found significantly different 
value of this index —1.80 suggesting that there may not be 
enough energy in microflares or nanoflares to heat the en- 
tire quiet corona. These last authors casted doubt on the 
hypothesis that the heating takes place in small scales, 
which according to them remains controversial. 

It is also worth mentionning that Benz & Kruckcr 
( |1998| ) and Berghmans et al. (1998) have noticed that the 
heating events occur not only on the magnetic network 
boundaries but also in the cell interiors. They came to the 
conclusion that these events have a quite similar nature. 

The idea of small scale intermittent behavior of dis- 
sipation was elaborated in theoretical models by several 
authors. Using an MHD approach, Einaudi & Velli (1982) 
investigated the dependence of the dissipated energy upon 
the scales. Different kinds of "fragmentation of dissipa- 
tion" were discussed by Vlahos (1982). Einaudi et al. 
Q1996D and Georgoulis et al. Ql99s| ) simulated 2D MHD 
systems driven by large scale forces and with diffusive dis- 
sipation. Their 2D simulations, driven by a pair of large 
scale vortices with random phases and amplitudes, showed 
that the spatial average of the dissipated power had a 
non-Gaussian statistics. After subtraction of the Gaussian 
component, the dissipated events were shown to have a 
scaling law distribution. 

The observed power-law distributions of the energy 
flux for the largest dissipative events, flares, and mi- 
croflares inspired the development of models based on the 
idea of Self-Organized Criticality (SOC) (Lu & Hamilton, 
19911 ; Lu et al., |199$ Vlahos et al . |1995| ; Georgoulis et al 



1995| ; Georgoulis & Vl ahos , Jl996| , |l998j j. Lu et al. ( [1993D 



and Georgoulis et al. (2001) compared the predictions of 
the models with the observations of the flares. SOC has 
appeared as a paradigm for slowly driven complex sys- 
tems, which exhibit power-law-type relaxation events and 



correlations of infinite range (Jensen, 1998). These mod- 
els are s imilar to the original "sandpile" model (Bak et 
al., 198S ) , but the sand is replaced by the magnetic field 



flux. The dissipation occurs through small scale reconnec- 
tion, when a "current," which is defined in terms of the 
difference between the field in a particular cell and the av- 
erage over all nearest neighbors, exceeds some threshold. 
In these models, a perturbation in a single cell can trig- 
ger quite a large cascade of reconnections. This results in 
a power-law energy distribution of dissipative events and 
power-law spatial correlations. However, one should take 
care interpreting these results because a small grid size 
can result in an artificial form of the PDF's obtained. A 



heuristic justification of such models based o n 2D MHD 
equations was proposed by Vassiliadis et al. (1995 ) . It is 
also worth mentioning tha t, usin g a shell model of MHD 
turbulence, Boffetta et al. ( 1999 ) have demonstrated that 
there can exist alternative reasons for a scaling law to ap- 
pear. 

Our model is also based on the idea of cellular au- 
tomata. Its difference with the previous ones consists in 
the assumption that not only the dissipation takes place 
in small scales but the magnetic field sources are also of 
small scale. This is in agreement with the idea of Benz 



and Kruckcr (1998) that the heating can take place on 
the level of chromosphere, implying that the magnetic 
field structures have the scales of the order or less than 
its height. This fe ature results from the analysis of Priest 
with co-authors ( 1998j ) who have shown that the heat- 



ing is quasi-homogeneous over the height of the magnetic 
loop. Other indication about the possible role of the small- 
scale sources follows from observations of Abramenko et 



al. ( 1999 ) of the evolution of the vertical component of the 
current helicity of an active region magnetic field during 
a transition from a low flaring state to an enhanced one. 
They found that the reorganization of the vortex structure 
of the photosphcric magnetic field occurs, when small scale 
vortices re-organize into larger scale vortices, suggesting 
that an inverse helicity cascade develops. 

Although the footpoints of the loops are known to 
move randomly, the statistical properties of their motions 
are not known in detail. These sources create strong cur- 
rents resulting in reconnection and flare-like events. Such 
eruptive events are related to magnetic energy releases, 
sudden changes of the magnetic field topology, acceler- 
ation of particles and plasma heating. Another possible 
source of energy dissipation and coronal heating is pro- 
vided by an anomalous resistivity resulting from current 
driven instabilities which are developed when the currents 
exceed a certain threshold value. It provides Joule- like dis- 
sipation in collisionless plasma. This results in a relatively 
slow heating, but does not give rise to particle accelera- 
tion as for the previous dissipation mechanism. Moreover 
it may occur everywhere in the turbulent shear flows and 
not necessarily on the boundaries of the magnetic field 
network. 

Another important feature that makes the difference 
from previous studies of CA type models is the homo- 
geneous distribution of the magnetic field sources. The 
problem to be addressed in the present work is whether 
the small scale magnetic fields generated in (or convected 
from) the photosphere or chromosphere (mainly at the 
boundaries between granulas, but also inside them) can 
make a significant contribution to the coronal heating. 

We are mostly interested in the scales typical for 
nanoflares and even in smaller ones. In this case, the small 
scale sources of the magnetic field and the energy dissi- 
pation can have comparable characteristic spatial scales. 
Thus the inertial interval as well as the energy cascade are 
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not so important in such a model as for the conventional 
Alfvenic or MHD turbulence. But the system considered 
can possess the property of "inverse cascade" , thereby in- 
fluencing the structure of the larger scale magnetic fields 
as is supposed to be the case in the SOC-type systems. In 
this paper we concentrate on the study of the dissipated 
energy flow, in particular, on the dependence of statistical 
characteristics of dissipated energy flow upon the charac- 
teristics of the magnetic field source, and on the physical 
mechanism of the magnetic field dissipation. 

2. Small scale driving and dissipation 




Fig. 1. A schematic representation of the procedure used 
in our model to simplify magnetic field configuration. 

We consider hereafter the sand-pile-type cellular au- 
tomata, where sand piles are replaced by a magnetic field 
flux. In the model, the relaxation process consists in the 
dissipation of the magnetic field energy by means of the 
dissipation of currents, when their magnitudes exceed a 
certain threshold value, or due to the reconnection pro- 
cess. 

In MHD approximation the evolution of the magnetic 
field is governed by 
OB 

= V x (u x B) + dissipative term, (1) 

where B is the magnetic field, it is the fluid velocity of the 
plasma. Without the dissipative term, this equation de- 
scribes the motion of the magnetic field lines frozen into 
the plasma. In this case, the reconnection is forbidden. 
When a finite conductivity a is taken into account in the 
form of Ohm's law j = a (E + u x B), where j is the cur- 



rent density, it gives rise to the diffusion effect. In Eq. (|lj) 
the effect is described by a dissipative term r]AB, where 
r] = 1/ er/io and /xo is the magnetic permeability of the free 
space. The diffusion time t v — I 2 /r) is generally large for 
all scales but small where it is comparable with the typical 
time for fast reconnection. In this case it has the same or- 
der of magnitude as an Alfven time t& — I/Va, where Va 
is the Alfven velocity. It is worth mentioning that at small 
scales other dissipative effects can become important (see, 
e.g., Biskamp 1993 ). 

In the present study, we use a simplified cellular au- 
tomata model mimicking Eq. (|]J). 

Assume that the curvature of the magnetic field lines 
is negligible. Then, in each cubic column only the upper 
and lower faces should be considered, i.e., magnetic field 
lines are straight and perpendicular to the surface. In this 
case the equation V • B = is automatically satisfied (see 
Fig. 0). As a result, the fields on the lower and upper faces 
are equal, thus only one face will be considered. The 2D 
surface is then split into cells. In most of the simulations 
presented in the paper, grids of 200 x 200 cells are used 
(unless otherwise stated). The boundary conditions are 
chosen to be periodic. 

As we have mentioned already in the introduction 
Krucker & Benz ( 1998) ) have found that the main heat- 
ing occurs in the small scale bright points. Making use 
of the multi- wavelength analysis Benz & Krucker (1999) 
have shown that the temporal sequence of observations of 
different wave emissions is similar to the one in the large 
scale magnetic loops. They came to the conclusion that 
the physical mechanisms of the energy release are similar 
in these two cases. It follows then that the dissipation pro- 
cess in the second case is also associated with the magnetic 
loops, but of a small scale. Moreover, they have noticed 
that the heating events occur not only on the boundaries 
of the magnetic network but in the interiors of the cells 
also. The comparative analysis of the model predictions 
for the plasma heating in the magnetic loop due to the 
distributed energy sourc e with observations performed by 
Priest with co-authors ( 1998 ) led to the conclusion that 
the heating is quasi-homogeneous along the magnetic loop. 
This means that the heating process does not occur in the 
close vicinity of the foot points but rather in the whole 
arc volume. If one will put together these facts, it follows 
that the characteristic spatial scale of the magnetic field 
loops which supply the magnetic field dissipated is of the 
same order as the characteristic scale of the dissipation. 
Thus we may conclude that not only the dissipative pro- 
cess, but also the sources have small characteristic length. 
Another conclusion from observations is that the sources 
are distributed quite homogeneously in space. This reason- 
ing leads us to the choice of the sources and dissipation 
mechanisms used in our model that we describe further. 

Three kinds of sources with slightly different statistics 
are considered. 
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2.1. Source terms 

We investigate the statistical behavior of the system 
driven by a random and turbulent unipolar/dipolar 
sources. These sources are used to model the effects of 
turbulent magnetic field convection described by the first 
term on the right-hand-side in Eq. ([!]). 

An unipolar random source. The simplest source 
of magnetic energy is an uncorrelated process of zero 
mean, the values (SB) = of which are randomly 
chosen from the set {—1,0,1}, all the values being 
cquiprobable. In each time step the action of the source 
consists in adding random numbers from the set men- 
tioned above to the previous values in the cell. The 
numbers are independently chosen for each cell. This 
procedure automatically ensures that (SB) = for 
each cell. 

— A dipolar random source. Such a source can be 
made dipolar by dividing the grid into two parts. For 
the positive and negative parts of the grid the ran- 
dom numbers (SB) = are chosen from the sets 
{-0.5, 0.5, 1.5} and {-1.5, -0.5, 0.5}, respectively. 

— A chaotic source. Turbulence is certainly not a com- 
pletely random process, and some of its aspects can 
be simulated using deterministic models. The Ulam 
map provides one of the simplest examples of a generic 
chaotic system with quadratic non-linearity (see, e.g., 
Frisch, 1995). The source SB in each cell evolves ac- 
cording to 

SB n+1 = l-2(SB n ) 2 , 

where the initial values of SBq are randomly chosen 
from the interval [—1,1]. In this case, all other 8B n 
belong to the same interval. The action of the source 
is similar to the above described for random sources. 

— A Geisel map source. In physical systems, the 
source term may depend on the value of B itself. When 
the dissipation is absent, the magnetic field in each cell 
evolves according to the map 



B„ 



f(Bn 



The initial values of Bq are randomly chosen from the 
interval [—0.5,0.5]. The systems like that are usually 



called Coupled Map Lattices (CML) (Kaneko, |1992[ ) 
rather than cellular automata. We use the source based 



on the Geisel map (Geisel & Thomae, 1984 ) shown in 
Fig. ||. The fixed points of this map that are defined 
by B n = f(B n ) are metastable. As a result, an inter- 
mittency develops in the system, i.e., the dynamics is 
slow and regular in the vicinity of the fixed points of 
the map and fast and chaotic otherwise. As a conse- 
quence, the map that describes the time evolution of 
the magnetic field in each particular cell exhibits the 
behavior similar to anomalous (non-Brownian) diffu- 
sion ,i.e. 

(B 2 )<xt a , a<l. 



The hypothesis has been suggested that when the tur- 
bulence consists of convective cells, magnetic field lines 
time evolution in each cell is similar to a subdiffusivc 
behavior (a < 1), which is, however, more complex 
than described above. 



c 
m 



-2 - 



-4 




Fig. 2. Graphical representation of the Geisel map (solid 
line). The fixed points of the map correspond to the in- 
tersections of the graph with the straight line B n+ i = B n 
(dashed line). 



2.2. Dissipation criteria 

The magnetic field dissipation provides the conversion 
of the magnetic energy into thermal energy and ensures 
the coupling between the magnetic field elements in our 
model. Phenomenologically, reconnection can be treated 
as a dissipation of small scale current sheets when the 
current density exceeds a certain threshold value (Somov 
& Syrovatsky, |l977t Syrovatsky, |l98l|, |1982|). 



If we neglect the displacement current, the current 
density can be calculated from Maxwell- Ampere's law 
VxB = jtioj, the finite-difference form of which can be 
written as 
Jx\ = J_( B(x,y)-B(x,y + S) 
jy J Sno\B(x + S,y)~B (x, y) 
where S is the grid increment. For simplicity we let S = 1 
and /^o — 1- Currents are thus computed as local gradients, 
and supposed to be carried along the borders between the 
cells. It is seen that the discrete analog of the current con- 
tinuity equation V • j = holds, i.e., the sum of incoming 
and outgoing currents is equal to zero at each node of the 
grid. 

The mechanisms for the current dissipation are of two 
types. 

— Anomalous resistivity, arises as the result of an in- 
stability when the electric current exceeds some critical 
value (which may be sometimes very small), and re- 
sults into an increase of the plasma's resistivity. This 
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phenomenon describes the phenomenon of pulse and 
energy exchange between electrons and ions, or be- 
tween different groups of particles of the same nature 
by means of plasma turbulence in collisionless plasma 
(Galeev & Sagdeev, |1979| , |1984| ). 

When the electric current in plasma exceeds certain 
threshold value plasma instabilities can be excited. 
As a result of the instability development the waves 
are generated (e.g. drift waves, or ion-acoustic, or low- 
hybrid waves). The development of the instability gives 
rise to the absorption by generated waves of a part of 
the electron's energy and pulse which is partly redis- 
tributed to ions. This results to a dissipation of the 
current quite similarly to Joule heating (see Appendix 
A for more details). Finally, it is worth noting that 
anomalous resistivity does not require any particular 
configuration of the magnetic field, besides the strong 
current, and may occur even in the presence of paral- 
lel magnetic fields pointing in the same direction, as 
for example in coronal holes or cells interiors. In our 
model we assume that the currents are completely an- 
nihilated whenever they exceed a certain threshold, 

\j\ > Jmax. 

Reconnection is generally understood as a relatively 
sudden change from one equilibrium state to another, 
implying a change in magnetic field's topology, accom- 
panied by a transformation of magnetic energy into 
energy of particules (Priest & Forbes, 20001 ). In its 
"primary form" the stationary reconnection process 
represents the dissipation of the magnetic field in the 
vicinity of the so called X-points of the magnetic field 
configuration or in the vicinity of a current sheet. That 
can be interpreted as the rising and "burnout" of thin 
small scale current sheets separating domains with op- 
positely directed magnetic fields. To mimic this process 
we assume that the following two conditions should be 
satisfied simultaneously, 



\j\ = \B-B'\ > j m 
B B' < 0, 



(2) 



where B and B 1 denote the magnitudes of the mag- 
netic field in neighboring cells. The conditions like that 
are supposed to hold for conventional plane reconnec- 



tion configuration with an X-point (Petchek, 1964; Sy 



rovatsky, 1981). Equation (0) results in the existence of 



currents that can significantly exceed the critical value 
(more details can be found in Appendix B). 
This is taken into account by requiring that magnetic 
fields in adjacent cells have opposite directions, which 
supposes the existence of a magnetic null point in be- 
tween, and corresponds to an unstable equilibrium fa- 
vorable to reconnection. 

This seems to be enough to distinguish between the 
two dissipation mechanisms in the framework of the cel- 



lular automata model imposing that the reconnection re- 
quires a special configuration (an X-point for instance), 
whereas the only criteria for anomalous resistivity is a 
current greater than a given threshold. Indeed, mag- 
netic reconnections corresponds to an important change 
in magnetic field's topology, whereas anomalous resistiv- 
ity does not. In real physical conditions of the corona both 
processes can be present. The threshold current for the 
anomalous resistivity is in general supposed to be much 
larger than for the reconnection under similar plasma pa- 
rameters. However our task in this work is confined in the 
study of the statistical properties of the energy dissipation 
dependence upon parameters. Closer comparison with real 
parameters in the corona will be the object of the more 
detailed study in the future and comes beyond the scope 
of this paper. 

The rules that we use to describe the dissipation pro- 
cess are based on the magnetic field dissipation, that 
means for us transformation to the particles energy or 
heating. For the sake of simplicity and pure formulation 
of the problem we consider here that all the magnetic field 
energy is transformed completely to the heating. The rules 
are as follows: when the current is annihilated, magnetic 
field values in both neighboring cells, B and B' , are re- 
placed by 1/2(5 + £>'), thus, the density of magnetic en- 
ergy dissipated in a single event is given by 



AE = - (B-B'f 
2 V ' 



1 .2 > 1 ,2 

2 J ~ 2 -'max" 



The procedure modeling the dissipation of currents is the 
same for both anomalous resistivity and reconnection. For 
each time step, the currents satisfying the dissipation cri- 
terion are dissipated till all the currents become subcrit- 
ical (or have the same sign in the case of reconnection). 
Then, we proceed to the next time step and switch on the 
source. Indeed, dissipative processes are supposed to be 
faster than the driving ones. The total dissipated energy 
is calculated as a sum over all the dissipated currents for 
the time step considered. 

To compare the effect of a single act of the magnetic 
field driver with that of the energy dissipation, one can 
analyze the spatial Fourier transform of the energy dissi- 
pation and the energy influx, 



16 
Z2T2 

^x^y 



-^diss(fc) 



16j 



2 

max 



V2 VI 
^x^y 



2 / ky 



.2 / ky 



where k Xl k y G [— 7T, tt]. The last equation holds for the dis- 
sipation of the y-component of current. For ^-component 
k x and k y should be interchanged. 

It is seen that the dependence of these two spectra on 
k x are quite different. For each fixed k y , the energy dissipa- 
tion vanishes at k x — and increases monotonously with 
the growth of |fc x |, while the source of the magnetic field 
is maximum at k x = and decreases with the increase in 
\k x \. It is worth noting that for large thresholds one act 
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of dissipation takes place after a large number of actions 
of the magnetic held inhux having on each step different 
phases, which are randomly distributed. This difference in 
the action of the magnetic held source versus dissipation 
introduces some intermediate scale in the fc-space where 
the action of the source is approximately compensated 
by the action of energy dissipation. This characteristic 
value of k is closely related to the characteristic correla- 
tion length of the magnetic held spatial distribution. The 
growth rate of the magnetic held energy and the effec- 
tive damping rate are non- vanishing almost everywhere in 
the fc-space, but the growth rate dominates for small fc's 
(large spatial scales) while dissipation dominates for larger 
k's (smaller spatial scales). Such a situation corresponds 
in terms of energy cascade to the "normal" direction of 
the energy flux, i.e., from large scales to smaller ones. It is 
also worth mentioning that using uniform driving results 
in an important difference between our model and conven- 
tional SOC models, where the extreme tenuousity of the 



driving is essential (Sornette et al., 1995). This tenuousity 
makes the driver nonlocal, in the sense that it depends on 
the state of the w hole system, as discussed by Vespigniani 
& Zapperi, |l998. 



is supposed to be of the order of 1, we have L ~ 300 m. 
This scale is significantly smaller than the spatial resolu- 
tion of modern experimental devices. Moreover, using the 
angular scattering measurements of the electron density 
fluctuations, the smallest scale that can be resolved in th e 
slow solar wind at 8i?© is 6 km (Woo & Habbal, 1997 ). 
Assuming that the linear structures expand radially as r, 
the structures at 1-R©, where dissipation is supposed to 
occur, are of the order of 1 km. 

Thus, until now, only macroscopic characteristics can 
be observed. The statistical microscopic models are aimed 
at reproducing the main features of these observations. In 
this paper we study the influence of the statistical prop- 
erties of the magnetic field source, of the type of dissipa- 
tion mechanisms and their thresholds on the macroscopic 
properties of the total flux of dissipated energy. This prob- 
lem is related to the other ones. In particular, can the lo- 
cal mechanisms of dissipation result in long-range spatial 
correlations? In self-organized systems, the appearance of 
such correlations give rise to power law distributions of 
dissipated energy. Such properties may be caused by spe- 
cific features of the source such as the deterministic chaos 
or non-Brownian diffusion. 



2.3. Characteristic spatial scales 

Dissipation mechanisms and their thresholds depend upon 
the parameters of the plasma of solar corona such as back- 
ground magnetic field, density, etc. Since our model is 
aimed at describing local regions in the corona rather 
than the corona as a whole, it is quite natural to as- 
sume that the same dissipation criterion can be applied 
for each cell of the grid. Thus the question about the 
characteristic sizes of the dissipated currents arises. The 
observations and theoretical studies show that the scale 
of current sheets can be smaller than 1 km. The smallest 
scales, of abou t 10 m, which are considered by Einaudi 
& Velli ( |l999|) , are associated with the current regions 
for the Petchek-type reconnection events. Let us estimate 
the characteristic scales of the dissipation events due to 
anomalous resistivity. Assuming that the resistivity is pro- 
vided by the ion sound instability that have a quite small 
threshold, we can easily obtain 

Ait 

|V x B\ ~ B/L > — n e ec s , 

where L is a characteristic width of the current sheet layer, 
c is the speed of light, B is the characteristic magnitude 
of the magnetic field, n e is the plasma density, and c s = 
(Te/rrii) 1 / 2 — (me/mi) 1 / 2 vr e is the ion-sound velocity, T e 
is the electron temperature, vx e is the electron thermal 
velocity, m e ^ are the electron and ion masses. Then 

AnnevTe \m e J Lu pl 
where u} p i is ion plasma frequency and (3 is the ratio of the 
kinetic to magnetic pressure. In the low corona, where (3 



3. Results 

Preliminary results concerned with the influence of the 
type of magnetic field energy dissipation on statistical 
properties of the total ra diatio n energy flux were presented 
by Podladchikova et al. ( 1999 ). It was shown that the dissi- 
pated energy has approximately normal distribution when 
the source of the magnetic field is random and the dissi- 
pation is provided by anomalous resistivity. This result 
was obtained under the condition that the current den- 
sity threshold for the dissipation to occur was moderate 
(imax = 5). Under the very same conditions, but for the 
dissipation due to reconnection, the non-Gaussian energy 
tails were observed as well as some other interesting fea- 
tures, in particular, large-scale spatial correlations of the 
magnetic field. This result seems to be quite natural be- 
cause in this case some currents may exceed the critical 
value but not dissipate when the magnetic fields have the 
same direction in the neighboring cells (see Eq. The 
currents may grow up to larger magnitudes and provide 
more intensive energy releases in a single dissipative event 
and sometimes longer chains of these events. The large en- 
ergy events are those that result in significant deviations 
from the Gaussian distribution. 



3.1. Transient and stationary states 

We have studied the dependence of the statistical prop- 
erties of our system upon the dissipation threshold j max 
ranging from 0.01 to 300. To perform statistical analysis 
correctly, the averaging procedures should be carried out 
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Fig. 3. Evolution of mean and variance of the magnetic 
field. For variance, observed are three stages — a linear in- 
crease, a nonlinear stage, and a stationary state. The re- 
sults are obtained for unipolar random source and anoma- 
lous resistivity dissipation with a threshold j max = 30. 

for the stationary state of the system. The transient time 
depends on the threshold. 

If the dissipation is absent, the field in each particular 
cell would follow a Brownian motion under the influence of 
the random source. As a result, the magnetic field in each 
particular cell, as well as the average over the whole box, 
would have a Gaussian distribution with a growing vari- 
ance, (B 2 ) oc t. The currents exhibit the similar behavior, 
the current variance growing as 

(j 2 )(t)*2{5B 2 )t, 

where for the random source (SB 2 ) =2/3. On the average, 
the current dissipates for the first time at a moment 



ts 



Jn 



2 < 5B 2 > ■ 



This relationship also gives a characteristic time between 
two dissipation events that occur in the same cell in a 
stationary state. To obtain a reliable statistics, the simu- 
lation time should significantly exceed ts- Depending on 
imaxj the simulations performed have 10 5 -10 6 time steps. 

Dissipation finally saturates the growth of the variance 
(see Fig. ||)b, and a stationary state is reached. The aver- 
age over the whole grid magnetic field B undergoes strong 
fluctuations but its time average is zero Fig. ||)a. The av- 
erage number n of dissipated currents at each time step 
can be estimated from the energy balance considerations. 
For a single time step, the energy input on a N x N grid 
is 

6E hl ~ N 2 (5B 2 ), 
while the dissipated energy is 



SE d 



In the equilibrium state 8E m ~ SE^iss, hence 
n ^ 4(SB 2 ) ^ 2 
N 2 " " 



J max 



ts 



(3) 



Assuming that these currents are uniformly distributed 
over the grid, the characteristic distance between them is 



I 



3.2. Grid size effects 



N 



1 H 



30 x 30 
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Fig. 4. Probability distribution function of the dissipated 
energy versus energy released (solid lines). Dotted lines 
represent the best fits by Gaussian distributions. The re- 
sults are obtained for unipolar random source and anoma- 
lous resistivity dissipation with a threshold j max = 5. The 
grid sizes are the following: a) 30 x 30; b) 50 x 50; c) 
64 x 64; d) 100 x 100. 

To investigate the grid size effects, we take grids con- 
sisting of 30 x 30-400 x 400 cells, the conditions of anoma- 
lous resistivity dissipation, and j max = 5. In this case we 
obtain the distributions that are very close to Gaussian for 
the largest grid. All the simulations were performed during 
1 6 time steps. The results are shown in Fig. |] for differ- 
ent sizes of the grid. For the grid 30 x 30, the distribution 
is nicely fit by a power law of index —3.1 for almost 3 en- 
ergy decades. As the grid size is increased, the bulk of the 
distribution becomes closer to a Gaussian one, while the 
high-energy tail retains a power-law shape. For the grid 
50 x 50 the power law shape of the tail is observed only for 
one energy decade, the estimate of the index being —2.9. 
For the grid 64 x 64 the power law tail is again shorter 
and the index is approximately equal to —4.3. For the grid 
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100 x 100 the distribution is practically indistinguishable 
from Gaussian in the whole energy range. 

Eventually, for large grid size the whole distribution 
becomes Gaussian. This means that the character of spa- 
tial correlations is changed. For small grids, spatial corre- 
lations decay rather slowly so that they extend over the 
whole grid. With the growth of the grid size, the expo- 
nential tails of the spatial correlations appear yielding the 
correlation length smaller than the grid size. The fact that 
the power-law distributions transform into Gaussian with 
the increase in the size of the system considered has al- 
ready been observed in some real sandpile experiments 
(Held et al. 



199C ) and in some forest-fire models (Grass- o- 



berger, |199l[ ). 

Now we proceed to the discussion of simulation results 
for various unipolar and dipolar sources and for both dis- 
sipation criteria using sufficiently large grids to avoid un- 
physical effects. 

3.3. Random unipolar sources 

Simulations on a 200 x 200 grid were performed for 2 x 10 5 
time steps and both types of dissipation (anomalous resis- 
tivity and reconnection) . The distributions of dissipated 
energy due to anomalous resistivity with the thresholds 
Jmax = 5, 30, 100 are shown in Fig. |[ Figure ^| represents 
the same PDF's but for dissipation provided by reconnec- 
tion. 

The reason for PDF shown in Fig. ||c to be multi- 
extremum is easy to understand. Since the dissipated 
energy in each event is of the form (_7 max + e) 2 , where 
< e < 2, for large j max , the peaks separated by j max a P~ 
pear. Assuming that e is a uniformly distributed random 
variable, we can get that the width of the peak number K 
is of the order of K(e 2 ). Hence, the discrete character of 
the PDF appears for the smallest values of the dissipated 
energy and large j max (see Fig. §b-c). 

From Figs. || and ^, one can see that the smaller cur- 
rent threshold (as compared to the source amplitude) re- 
sults in a PDF of dissipated energy that is close to nor- 
mal. With the growth of j max , a high energy nonthermal 
tail appears. Although the average value of the dissipated 
energy does not depend upon j max , the deviations from 
Gaussian distribution in the tail increase with j max and 
are more pronounced for dissipation due to reconnection. 
For the dissipation due to reconnection, all three distribu- 
tions have non-Gaussian tails that can be approximated 
by power-law. With the increase of j max the tails become 
more extended and the index of the power-law distribu- 
tion decreases. This signifies that the stronger deviations 
from Gaussian distribution appear for larger values of the 
threshold. 

The formation of such a tail is not related to the in- 
crease of the correlation length, because the largest mag- 
netic field correlation length is about 20 for moderate 
threshold, j max = 5, in the case of reconnection-like dis- 
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Fig. 5. Probability distribution function of the dissipated 
energy versus energy released (solid lines). Dotted lines 
represent the best fits by Gaussian distributions. The re- 
sults are obtained for unipolar random source and anoma- 
lous resistivity dissipation. The threshold currents j max 
arc the following: a) j max = 5; b) j max = 30; c) j max = 
100. 



sipation. Let us also notice that the average correlation 
length is smaller when anomalous resistivity dissipation is 
consid ered a nd decreases as j max increases (Podladchikova 
1999)) . 



et al., 

The average number of dissipated currents decreases 
with the increase of the threshold. Indeed, it easily seen 
from Eq. (||) that for j max = 1 we have 



n 

iV2 



4 
3' 



i.e., each current is dissipated at each time step and dissi- 
pative events occur almost everywhere on the grid thereby 
creating long range correlations. For large j max , the ra- 
tio n/N 2 becomes small, e.g., for jmax = 100 we have 
n/N 2 ~ 3 • 10~ 4 that corresponds to the uncorrelated dis- 
sipative events with relatively small overlapping. 

Thus, for large j max the dissipation is much faster than 
the action of the magnetic field source, i.e., the conditions 
of time scale separation are better satisfied. 
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Fig. 6. Same as for Fig. but for the reconnection-like 
dissipation. 

3.4- Ulam map source 

Despite the deterministic nature of the magnetic field in- 
crement SB, without dissipation the magnetic field in each 
cell exhibits a Brownian motion with Gaussian statistics 
after a large number of steps. Indeed, each SB has the 
same PDF with a finite variance. Thus, the dependence 
of the PDF upon the critical current and the dissipation 
mechanism have manifested itself in similar tendencies as 
random sources do. The PDF's for moderate critical cur- 
rent, jmax = 5, are close to Gaussian, for j max = 30 the 
small deviations from Gaussian distribution begin to ap- 
pear for high energies, and for j max = 100 there is clear 
evidence of the presence of the high energy suprathermal 
tail. The tail is much better pronounced for the dissipation 
due to reconnection (Figure 0). 

3.5. Geisel map lattice 

Another source of the magnetic field is provided by the 
Geisel map. It was used as a source of the field rather 
than the increment of the field. As it was discussed above, 
this map exhibits the time evolution in each cell quite sim- 
ilar to anomalous (non-Brownian) diffusion. The growth 




4 



2x10 4 4x10 4 6x10 4 



E 

Fig. 7. Probability distribution function of the dissipated 
energy versus energy released (solid lines). Dotted lines 
represent the best fits by Gaussian distributions. The re- 
sults are obtained for Ulam map source and j ma x = 30. . 
The dissipation laws are the following: a) anomalous re- 
sistivity dissipation; b) reconnection-type disispation. 
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Fig. 8. Probability distribution function of the dissipated 
energy versus energy released (solid lines). Dotted lines 
represent the best fits by Gaussian distributions. The re- 
sults are obtained for Geisel map source and anomalous 
resistivity dissipation. Threshold current j max = 30. 



of the magnetic field variance is slower than for Brown- 
ian motion. Such a behavior of the walker would be called 
subdiffusive. As a result, the variance is smaller than for 
the random source. This causes some differences observed 
in the behavior of the PDF of the dissipated energy. If 
the dissipation is provided by anomalous resistivity, the 
PDF's obtained with this map and the random source are 
quite similar except the width of the distribution. Fig- 
ure H represents the PDF obtained for Geisel map and 

Jmax — 30. 

The smaller current threshold (as compared to the 
source amplitude) results in a PDF of dissipated energy 
that is close to normal. With the growth of j max , a high en- 
ergy nonthermal tail appears. The deviations from Gaus- 
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Fig. 9. Probability distribution function of the dissipated 
energy versus energy released (solid lines). Dotted lines 
represent the best fits by Gaussian distributions. The re- 
sults are obtained for Geisel map source and dissipation 
due to reconnection. Threshold current is j max = 30. 



sian distribution in the tail increase with j max . With the 
increase of j ma x the tails become more extended. 

In this case the range of energies is a slightly less than 
for the Ulam map, the distribution being quite jagged, all 
the other features of the distributions seem to be similar. 

Figure || shows the PDF obtained for the same Geisel 
map source, dissipation due to reconnection, and j max = 
30. A pronounced tail of the distribution is evidenced 
thereby confirming the tendencies already observed for 
other types of sources. A difficulty that one encounters 
when analyzing such PDF's is illustrated in Fig. |[ In 
Fig. ||a, the PDF is represented in semi-logarithmic scale 
and the distribution seems to have an exponential tail. In 
Fig. ^|b the same distribution is shown in the log-log repre- 
sentation, where the same tail can be treated as a power- 
law- type distribution with the index —3.1. This problem 
is related to the fact that the tail extends over a rather 
short range of about one energy decade. This problem will 
be tr eated in more detail elsewhere (Podladchikova et al. 
2001) . 




Fig. 10. Magnetic field structure observed in the simula- 
tions with the dipolar random source, anomalous resistiv- 
ity dissipation, and a threshold current j max = 5. 

localized current layer in the transition region between the 
two parts of the system. The grid size is 100 x 100. 

The dissipation provided by reconnection is not always 
able to stabilize the field growth. However, for dissipation 
due to anomalous resistivity a stationary state always ex- 
ists. Further we consider only the case with anomalous 
resistivity dissipation. 

The energy of the dipolar structure grows faster than 
in the unipolar case, and a stationary state is quickly 
reached. The dipolar structure of the magnetic field is 
presented in Fig. [n| This structure is stationary. There- 
fore the spatial correlations of magnetic field are long and 
quasi-stationary. 

Such a system differs from the unipolar one with the 
random source. We observe, as could be expected, that the 
dissipation events are more intense within the transition 
region where the magnetic field changes its sign. In addi- 
tion, in the dipolar case the position of the PDF peak value 
on the energy axis is significantly shifted towards high en- 
ergies with the increase of the threshold. Such a property 
can be explained by the presence of the localized current 
layer in the vicinity of the neutral sheet. The characteris- 
tic dependence of the PDF's upon the dissipated energy 
in dipolar case displays a supra-thermal high energy tails 
which were not observed in the unipolar case with anoma- 
lous resistivity (see Fig. |ll|). The characteristic indices for 
jmax = 50, 100, and 230 are approximately equal to —1.7, 
— 1.9, and —2.3, respectively. For j max — 230, the time 
series of the dissipated energy is shown in Fig. O. 



3.6. Dipolar source 

When the effects due to the bi-polarity are modeled, the 
grid is split into two equal parts. For one half of the grid, 
the value of the magnetic field increment is chosen ran- 
domly from the set {—1.5, —0.5,0.5}, for another half the 
set {1.5,0.5,-0.5} is used. The global inhomogeneity of 
the dipolar source results in the formation of a strongly 



4. Discussion and conclusions 

To study coronal heating due to dissipation of small-scale 
current layers, we have performed a statistical analysis of 
a simple cellular automata model. Its principal difference 
with previous ones is that the system is driven by small 
scale homogeneously distributed sources acting on the en- 
tire grid for each time step. The idea to consider small 
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Fig. 11. Probability distribution function of the dissi- 
pated energy versus the energy released (solid lines). Dot- 
ted lines represent the best fits by Gaussian distributions. 
The results are obtained for dipolar random source and 
anomalous resistivity dissipation. The duration of calcu- 
lations was 2 x 10 5 time steps. The threshold currents 
jmax are the following: a) j max = 5; b) j max = 50; c) 
jmax = 100; d) j max = 230. 

scale sources is similar to the one proposed by Benz & 
Kruckcr that the heating occurs on the level of the chro- 
mosphere, thus the magnetic field structures, dissipation 
of which supplies the energy for the heating, are also of a 
small scale. 

The magnetic field sources we use are either of the 
following types: 

— random source, 

— deterministic chaotic map (Ulam map), both for the 
magnetic field increment. 
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200 400 600 800 
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Fig. 12. Fragment of time series of the dissipated energy, 
for jmax = 230, dipolar random source, anomalous resis- 
tivity dissipation. 

— Geisel map (coupled map lattice) that describes the 
time evolution of the magnetic field in each cell similar 
to anomalous (non-Brownian) diffusion. 

We consider two mechanisms of small scale dissipation 
to occur: 

— the first is used to model anomalous resistivity dissi- 
pation, 

— the second is used to model the local reconnection. 

The first one is supposed to be similar to Joule dissipa- 
tion but in collisionlcss plasma. It relies on the criterion 
that any local current whose magnitude exceeds a pre- 
determined threshold value must dissipate. 

In the second case we assume that the dissipation oc- 
curs when the above mentioned condition is satisfied (this 
docs not mean that the threshold in real physical system 
should be the same in two cases) and, in addition, the 
magnetic fields vectors in neighboring cells should have 
opposite signs. This mimics the presence of the null point 
of the magnetic field in the reconnection-type configura- 
tions. 

The characteristic under investigation is the total en- 
ergy of all simultaneously dissipated currents. 

Our observations can be summarized as follows. 

Small scale magnetic field sources and localized energy 
dissipation mechanisms can result in large scale correla- 
tions of the magnetic field. In our calculations, the energy 
sources dominate in larger scales while the dissipation in 
the smaller ones. Thus, the system behaves as having the 
energy cascade from large scales to smaller ones. However, 
during the evolution of the system, large scale correlations 
are formed with a characteristic length significantly larger 
than that of the source and dissipation. 

For all three types of the magnetic field sources consid- 
ered, the obtained PDF of the dissipated energy is close 
to Gaussian distribution when relatively small threshold 
currents, jmax < 10, are chosen. The dependence of the 
PDF of the dissipated energy upon statistical properties 
of the source for all three types of the sources considered is 
rather weak and requires an additional thorough analysis 
to be performed. 
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In the case of reconnection type dissipation, the de- 
viations from Gaussian distribution are stronger than 
for anomalous resistivity dissipation. For large values of 
the threshold current density, we observe a high-energy 
suprathermal tail which have the shape similar to power- 
law distribution. The time-averaged spatial correlations 
are exponentially decaying. In this sense our model can 
not be considered as a self-organized critical system. How- 
ever, sometimes long range (power-law type) correlations 
are observed. 

To model a dipolar global magnetic field structure 
we used inhomogeneous distribution of the magnetic field 
source. It has small average positive value in the one part 
of the grid and negative in another with the summary 
average equal to zero. The system does not reach any sta- 
tionary state in the case of the reconnection type dissipa- 
tion. In the case of anomalous resitivity it saturates and 
the stationary state is established. In this last case, the de- 
viations from the Gaussianity become stronger with the 
increase of the threshold. This effect is associated with 
the decrease of the thickness of the current layer where 
the dissipation is concentrated. 

The statistical analysis of frequency distribution of 
such heating events shows always larger absolute value 
of the exponent for the more energetic reconnection pro- 
duced events than for the anomalous resistivity heating 
events. This tendency seems to be similar to the one found 
by Benz & Krucker, 199§| who studied the emission mea- 
sure increases. They pointed out that the statistical prop- 
erties of the faint events that occur in the intra-cell re- 
gions of the quiet Sun manifest quite small deviations from 
Gaussian distributions, while supposed nanoflares that are 
associated with the network boundaries have more promi- 
nent enhancements and stronger deviations. Taking this 
into account we come to the conclusion that the quanti- 
tative difference between faint and strong heating events 
reported in (Krucker & Benz, 200C ) can probably be ex- 
plained by two different mechanisms of magnetic field dis- 
sipation that we used in our model. 

However, our current work represents only the first 
step in the development of a model. In particular, many 
elements of its should be better adapted to real physi- 
cal processes. Further effects should be included for direct 
comparison with the observations, such as more detailed 
description of respective time scales of the instabilities 
or particle acceleration, thermal/non-thermal processes or 
outgoing flows escaping from the reconnection region. For 
instance one can assume that in the case of the anomalous 
resistivity the dissipation of the current can be partial in 
spite of the complete, the two dissipation mechanisms can 
be combined, anomalous resistivity can co-exist with the 
reconnection. One can more carefully take into account the 
energy re-distribution between different particle species as 
we have noted in the Appendices. Thus, on this stage we 
can only conclude that some tendencies of our model seem 



to be in a good qualitative agreement with the experimen- 
tal ones. 
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Appendix A: Anomalous resistivity 

In this appendix, we explain in more details the physics 
of the "anomalous resistivity" mentioned in our paper. 
The idea is based on the possibility of the exchange of 
the pulse and energy between the electrons and the ions 
due to the instability development and the appearance 
of the turbulent state in the plasma. The "conventional" 
expression for the conductivity reads 

a = (ne 2 ) / (m e v) , 

where n is the plasma density, e, m e - electron charge and 
mass respectively, v is the collision frequency of the elec- 
trons with ions. In collisionless plasma the collision fre- 
quency is negligible and the collisional conductivity is in- 
finite. However the electrons (who carry the major part 
of the electric current) can excite collective oscillations of 
the electrons as well as ions, and transfer part of their im- 
pulse and energy to these oscillations. This results to an 
anomalous loss of the electron pulse and energy, and, con- 
sequently, to the decrease of their directed velocity, i.e. the 
decrease of the current. Such a process can be character- 
ized by an "effective collision frequency" f e //. To calculate 
this characteristic frequency one can consider the pulse 
conservation law in the system that consists of electrons 
and waves. To this end one should evaluate this effect as 
the action of the friction force that decelerates the elec- 
tron flow. The pulse loss per unit time can be written as 
follows: 

f e //m e nUd = -F fr 

where Ud is the directed velocity of the electron flow that 
carries the current, and Ff r the friction force that acts on 
the electrons. The same decrease of the electron pulse can 
be estimated as the increase of the pulse of the waves emit- 
ted by the electrons taking into account that the change 
of the pulse of the waves is described by the following 
expression: 



dP w 
dt 



lkW k 



k d 3 k 

Uk (2tt) 3 



where 7^ is the linear increment of the instability and Wk 
is the spectral energy density of waves. The assumption 
that the pulse of the system consisting of electrons and 
waves is conserved results in the equality of these two 
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expressions, thus the effective collision frequency can be 
defined as: 



m e nUl 



(U d k) d 3 k 

u k (2tt) 3 



Thus the estimate of the level of the wave turbulence tak- 
ing into account the nonlinear saturation mechanism al- 
lows to solve the problem. 

When the current flows in the direction perpendicular 
to the magnetic field two major types of the instabilities 
that give rise to the anomalous resistivity are modified 
Buncman instability and the instability with respect to the 
generation of the so called electron-acoustic modes. The 
final effect of the anomalous turbulent resistivity is the 
transfer of the energy from the electrons to ions, because 
the damping of these waves takes place mainly due to their 
interaction with ions. 

In the particular case of a current perpendicular to 
the magnetic field, two major types of instabilities that 
result in the anomalous resistivity are modified Buneman 
instability and the instability which generates the so called 
electron-acoustic waves. 

An important physical characteristics of the anoma- 
lous resistivity phenomenon is the ratio of the energy dis- 
sipated by ions and to the energy dissipated by the elec- 
trons. If the characteristic drift velocity of the electrons 
with respect to ions is Vd, and we have the knowledge 
of the spectral characteristics of the wave spectrum ex- 
cited, let us assume, for instance that the characteristic 
frequency of the waves excited due to the instability is 
Uk, and that characteristic wave vector of the unstable 
waves is k, then this ratio can be estimated by: 

dE 1 dE ir /7kVF k ^d 3 k ( kVd ) 
dt ' dt 



J 7 kWkd 3 k 



for typical instabilities, where this ratio is approximately 
1. 

So, the final result of the anomalous resistivity is quite 
similar to the Joule heating of the ion component of the 
plasma. The energy dissipation that heats the plasma can 
be represented as: 

Q =f/°eff 

This phenomenon was observed e xperi mentally in the lab- 
oratory plasma (Eselevich et al., 1971). 



Appendix B: Reconnection 

There are several differences between reconnection process 
and anomalous resistivity. One of them is related with the 
change of the magnetic field topology in the first case while 
in the second there are only quantitative variations of the 
basic characteristics of the magnetic field configuration. 

Another difference is that while the diffusion processes 
in the first case provide heating (comparable for ions and 



electrons, see appendix), reconnection converts magnetic 
energy mostly into particle (mainly ions) acceleration. 
Then energetic beams may provide heat, but this is only 
an indirect consequence of the reconnection. One more but 
quite important for observations is a difference in time 
scales. The reconnection is supposed to be a rapid energy 
release while the anomalous resistivity is relatively slow 
diffusive process. 

Moreover, as it was shown analytically and in com- 
puter simulations, the reconnection process can under cer- 
tain conditions look like an explosive event. During the 
reconstruction of the magnetic field topology the compo- 
nent of the magnetic field perpendicular to the background 
magnetic field and to the direction of the current can grow 
so rapidly that its variation in time is quite similar to an 
explosion. The spatial temporal dynamics of the magnetic 
field can be described by the following expression: 



where 



L 



B z (t) 



B = Bq x tanh I — \ e x + B z (t) sin (kx) e z , 



Bq z 



1 t / T eX pl 



(see Galeev, 1984 for more details). Here Bq x is the mag- 
nitude of the background surrounding magnetic field, that 
is supplied by the current carried along the y axes, r exp i is 
the characteristic time of the magnetic and electric fields 
variations. It has opposite sign from two sides of the cur- 
rent layer. B z is the initial amplitude of the perturbation 
of the normal component of the magnetic field that is sup- 
posed to be periodic along the x axes with the wavelength 
A = 27rk). Its time dependence represents the growth of 
the amplitude up to infinity for the finite time. However 
this formal solution is valid only when the amplitude of 
this perturbation is smaller than the background field, 
thus the singularity is not so important. This explosive 
growth of nonlinear perturbations results in the similar 
increase of the inductive electric field. It does not mean 
that the electric and magnetic fields become infinite, just 
the nonlinear stage develops in such a rapid way till the 
normal component of the magnetic field becomes compa- 
rable with the Box- This process results in the rapid ac- 
celeration of electrons and ions by the inductive electric 
field in the region where the particles are unmagnetized. 

Another observational feature of the reconnection pro- 
cess consists in the presence of the macroscopic fluxes 
around the reconnection site with the characteristic ve- 
locities of the order of the Alfven speed estimated using 
background magnetic field in the vicinity of it. These flows 
can give rise to the replenishing of the lower density re- 
gions by the material from the reconnection site if the den- 
sity there is larger or even comparable with the density of 
the surrounding plasma. A part of the energy released is 
supposed to be transformed into the kinetic energy of ac- 
celerated charged particles. It is known that the energy 
transferred to particles growth with the increase of their 
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mass as M 1 / 3 , the heaviest one s bein g the most effectively 
accelerated (Vekstein & Priest, 1995 ). In their turn, these 
particles can generate the electromagnetic radiation pro- 
viding experimental signatures of the heating events. 
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